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We  describe  a  recently  discovered  stable  planar  surface  of  silicon,  Si(  1 14).  This  high-index  surface, 
oriented  19.5°  away  from  (001)  toward  (111),  undergoes  a  2  X  1  reconstruction.  We  propose  a 
complete  model  for  the  reconstructed  surface  based  on  scanning  tunneling  microscopy  images  and  first- 
principles  total-energy  calculations.  The  structure  and  stability  of  Si(l  14)-(2  X  1)  arises  from  a  balance 
between  surface  dangling  bond  reduction  and  surface  stress  relief,  and  provides  a  key  to  understanding 
the  morphology  of  a  family  of  surfaces  oriented  between  (001)  and  (114).  [S003 1-9007(96)007 18-1] 

PACS  numbers:  68.35.Bs,  61.16.Ch,  73.61. Cw 


The  (001)  and  (111)  surfaces  of  silicon  are  its  primary 
growth  and  cleavage  planes,  respectively.  Surfaces  with 
orientations  between  these  two  planes  are  both  intrinsically 
interesting  and  technologically  important.  For  these  inter¬ 
mediate  orientations,  surface  morphology  depends  strongly 
on  the  orientation  angle.  For  example,  small  misorien- 
tations  away  from  the  low-index  planes  produce  vicinal 
surfaces  consisting  of  terraces  separated  by  steps  [1-3], 
whereas  larger  angles  often  lead  to  sawtoothlike,  grooved 
surfaces  [4,5].  Between  (001)  and  (111),  only  two  orien¬ 
tations  are  known  to  form  planar  surfaces  with  stable  re¬ 
constructions:  Si(l  13)  [6,7]  and  Si(5  5  12)  [8].  Recently, 
(114)  surfaces  have  been  observed  both  on  cylindrical  Si 
samples  [9]  and  within  etch  pits  formed  on  Si(001)  [10]. 
In  addition,  sawtoothlike  structures  composed  of  (114)  and 

(113)  facets  have  been  observed  on  samples  oriented  be¬ 
tween  (114)  and  (113)  [5,9].  These  results  suggest  that 
Si(  1 14)  is  a  planar  surface  that  is  thermodynamically  stable 
against  faceting. 

The  (114)  plane  can  be  considered  a  highly  vicinal  (001) 
surface,  oriented  0  =  19.5°  from  (001)  toward  (1 1 1).  The 
structure  of  the  ideal  bulk-terminated  (114)  surface  is 
an  alternating  sequence  of  one-  and  two-unit-cell  wide 
(001)  terraces  separated  by  double-layer  steps,  as  shown 
in  Fig.  1(a).  This  basic  terrace-plus-step  morphology  also 
describes  real  surfaces  over  a  large  range  of  orientation 
angles:  from  slightly  vicinal  (001)  surfaces  with  large 
variable-width  terraces  to  steeper  orientations  consisting 
of  three-unit-cell  wide  terraces  separated  by  double-layer 
steps  ( 6  =  11.4°)  [11].  It  is  interesting  to  ask  at  what 
orientation  angle  this  morphology  breaks  down,  and  for 
what  reasons.  For  example,  the  Si(  113)  surface  ( 0  = 
25.2°)  has  a  stable  3X2  reconstruction  that  cannot  be 
simply  described  as  the  expected  sequence  of  one-unit-cell 
terraces  and  double-layer  steps  [7].  For  a  surface  with 

(114)  orientation,  one  is  therefore  led  to  consider  three 
possibilities:  (i)  a  periodic  surface  with  a  terrace-plus-step 
morphology;  (ii)  a  periodic  surface  without  a  terrace-plus- 
step  morphology,  similar  perhaps  to  Si(113)-(3  X  2);  or 
(iii)  a  nonperiodic  (e.g.,  grooved)  surface.  In  this  Letter  we 
use  scanning  tunneling  microscopy  (STM)  to  definitively 
establish  that  Si(114)-(2  X  1)  is  a  planar  surface  with 


a  stable  terrace-plus-step  reconstruction.  We  propose  a 
complete  structural  model  for  the  reconstructed  surface, 
and  use  first-principles  electronic-structure  methods  to 
provide  strong  theoretical  evidence  in  support  of  this 
model.  The  structure  of  Si(  1 14)  reveals  it  to  be  a  natural 
extension  of  vicinal  (001);  indeed,  it  marks  the  end  point 
of  a  family  of  surfaces  having  closely  related  structure. 

The  clean  Si(  114)  surface  [12]  appears  in  constant- 
current  STM  images  as  a  well-ordered  periodic  array 
of  row  structures  oriented  along  the  [110]  direction 
[Fig.  1(b)].  The  period  of  this  structure  in  the  [221] 
direction  is  approximately  16  A,  equal  to  the  length  of 
the  bulk-terminated  (114)  unit  cell.  Along  each  row,  the 
period  in  the  [110]  direction  is  7.7  A,  twice  that  of  the 
bulk-terminated  unit  cell.  Atomically  resolved  images 
reveal  three  distinct  types  of  row  structures,  labeled  I,  II, 
and  III  in  Fig.  1(c).  The  phase  relationship  along  [110] 
between  these  row  structures  varies  over  different  regions 
of  the  surface,  giving  rise  to  areas  of  local  2  X  1  as  well 
as  c(2  X  2)  periodicity  [see  Fig.  1(c)]. 

Our  proposed  structural  model  for  Si(114)-(2  X  1)  is 
closely  related  to  the  structure  of  vicinal  Si(001).  The 
Si(001)  surface  undergoes  a  reconstruction  in  which  the 
surface  atoms  pair  up  and  bond  to  reduce  the  number 
of  surface  dangling  bonds,  forming  parallel  rows  of 
dimers.  For  vicinal  surfaces  with  0  greater  than  4°- 
5°,  the  great  majority  of  resulting  (001)  terraces  consists 
of  dimers  oriented  parallel  to  the  step  edges  (/(-type 
terraces),  with  these  /(-type  terraces  separated  by  double¬ 
layer  steps  ( Db  steps)  [1].  Chadi  first  proposed  that 
Db  steps  are  rebonded;  i.e.,  they  incorporate  another 
row  of  atoms  at  the  step  edge  to  reduce  the  number  of 
dangling  bonds  [13].  This  basic  morphology — /(-type 
terraces  separated  by  rebonded  DB  steps — has  also  been 
observed  in  STM  studies  of  curved  Si  samples  having 
regions  with  local  orientation  6  ~  9°-ll°,  where  the 
(001)  terrace  width  is  only  three  or  four  unit  cells  [11]. 
We  propose  that  a  variant  of  this  basic  morphology  also 
describes  the  very  narrow  terraces  found  on  Si(l  14). 
Specifically,  our  structural  model  for  Si(114)-(2  X  1), 
shown  in  Fig.  2,  is  an  alternating  sequence  of  single  - 
and  double-width  /(-type  (001)  terraces,  separated  by 
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FIG.  1.  (a)  Side  view  of  the  Si  crystal  lattice,  showing 

one  unit  cell  of  the  bulk-terminated  (114)  surface.  Dashed 
lines  emphasize  the  alternating  single-  and  double- width  (001) 
terraces  separated  by  double-layer  (DL)  steps,  (b)  Filled-state 
STM  gray-scale  image  (210  X  180  A)  of  the  clean  Si(114) 
surface,  (c)  Atomic-resolution,  filled-state  image  (65  X  54  A) 
highlighting  the  three  types  of  observed  row  structures.  Unit 
cells  are  drawn  to  indicate  two  types  of  local  symmetry: 
2X1  (left)  and  c(2  X  2)  (right).  The  2X1  unit  cell  is 
16.3  X  7.7  A. 

rebonded  Db  steps  (stepping  up  from  the  double-width 
terrace)  and  nonrebonded  DB  steps  (stepping  down  from 
the  double-width  terrace).  The  presence  of  both  rebonded 
and  nonrebonded  steps  provides  an  important  mechanism 
for  stress  relief  within  the  unit  cell,  as  discussed  in  de¬ 
tail  below. 

This  structural  model  is  strongly  supported  by  the  re¬ 
sults  of  first-principles  calculations  of  equilibrium  surface 
geometry  and  surface  energy,  and  the  resulting  simulated 
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FIG.  2.  Top  and  side  views  of  proposed  model  for  Si(l  14)- 
(2  X  1)  reconstruction,  with  fully  relaxed  geometry.  Labeled 
structures  correspond  to  the  STM  features  seen  in  Figs.  1  and 
3.  A  dimer  plus  two  nonrebonded  step-edge  atoms  is  called 
a  “tetramer.” 

STM  images.  The  calculations  were  performed  in  a 
supercell  slab  geometry  (inversion  symmetric)  with  ten 
layers  of  Si  and  a  vacuum  layer  equivalent  to  six  Si  layers. 
Full  structural  relaxation  was  performed  using  the  Corn¬ 
ing  electronic-structure  code  of  Teter,  Payne,  and  Allan 
[14],  which  solves  the  Kohn-Sham  equations  in  the  local- 
density  approximation  (LDA)  with  a  plane-wave  basis 
and  norm-conserving  pseudopotentials  [15].  The  kinetic- 
energy  cutoff  was  10  Ry,  and  a  single  k  point  was  used 
throughout.  Structural  relaxation  was  performed  on  all 
the  atoms  until  the  rms  force  was  less  than  0.1  eV/A  [16]. 

The  fully  relaxed  surface  geometry  for  our  proposed 
model  includes  three  prominent  row  structures  per  unit 
cell,  corresponding  to  the  three  types  of  rows  seen  in  the 
STM  images.  Row  I  in  the  images  arises  from  the  line  of 
rebonding  atoms  at  the  rebonded  DB  step  edge;  row  II 
arises  from  the  row  of  tetramers  (i.e.,  the  combination 
of  the  dimers  on  the  double-width  terrace  and  the  nonre¬ 
bonded  Db  step-edge  atoms);  and  row  III  arises  from  the 
dimers  on  the  single-width  terrace.  The  correspondence 
between  theory  and  experiment  is  examined  in  detail  in 
Fig.  3.  Figures  3(a)  and  3(b)  show  STM  images  for  filled 
and  empty  states,  respectively,  from  the  same  surface  area; 
Figs.  3(c)  and  3(d)  show  numerically  simulated  STM  im¬ 
ages  (constant-current  surfaces)  calculated  by  integrating 
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FIG.  3.  (left)  STM  images  (33  X  31  A)  of  Si(l  14)  mea¬ 
sured  for  (a)  filled  states,  1.2  V;  (b)  empty  states,  —1.7  V. 
(right)  Simulated  STM  images  calculated  for  (c)  filled  states; 
(d)  empty  states.  Black  dots  indicate  projected  atom  posi¬ 
tions.  The  simulated  images  are  for  a  surface  with  c(2  X  2) 
symmetry. 

the  local-state  density  over  filled  and  empty  states,  respec¬ 
tively,  near  the  Fermi  level.  In  the  filled-state  STM  image 
[Fig.  3(a)],  rows  I,  II,  and  III  all  show  2ao  periodicity  along 
[110]  (no  =  3.84  A  is  the  surface  lattice  constant):  rows  I 
and  III  each  appear  as  an  array  of  maxima  with  2gq  spac¬ 
ing,  whereas  row  II  has  a  more  variable  appearance.  In  the 
empty-state  image  [Fig.  3(b)],  all  three  rows  also  exhibit 
2gq  periodicity,  but  the  features  are  more  diffuse  within 
each  row. 

The  general  appearance  of  the  rows  is  accurately  repro¬ 
duced  in  the  simulated  STM  images  based  on  our  model. 
In  the  simulated  filled-state  image  [Fig.  3(c)],  row  I  is  most 
prominent,  and  the  expected  2gq  periodicity  is  apparent 
in  all  three  rows,  in  good  agreement  with  the  STM  data. 
Along  row  I,  weak  maxima  are  visible  between  the  domi¬ 
nant  peaks  spaced  2ao  apart;  these  maxima  are  only  occa¬ 
sionally  resolved  in  STM  images.  Row  II  has  a  zigzaglike 
internal  structure,  consistent  with  the  variable  appearance 
observed  experimentally.  In  the  simulated  empty-state 
image  [Fig.  3(d)],  the  2gq  periodicity  is  still  evident  in 
row  II,  but  is  much  weaker  in  rows  I  and  III,  as  observed  in 
the  experimental  image.  Finally,  the  experimental  obser¬ 
vation  of  areas  with  local  2X1  and  c(2  X  2)  symmetry  is 
also  consistent  with  our  model.  We  have  calculated  total 
energies  of  both  reconstructions  and  find  differences  only 
at  the  level  of  1-2  meV/A2. 

Several  of  the  topographic  features  described  above  are 
a  consequence  of  a  structural  distortion  that  occurs  within 
each  row.  For  example,  at  equilibrium  the  dimers  that 
constitute  row  III  are  each  tilted  (buckled)  so  that  the  two 


dimer  atoms  have  a  theoretical  height  difference  of  Az  = 
0.17  A.  This  buckling  occurs  on  flat  Si(001)  as  well, 
and  is  energetically  favored  by  the  transfer  of  an  electron 
from  the  lower  to  upper  atom  of  the  dimer.  A  similar 
buckling  occurs  along  the  line  of  rebonding  atoms  that 
constitute  row  I  (Az  =  0.18  A).  Both  of  these  bucklings, 
although  small,  give  rise  to  a  large  corrugation  of  the 
constant-current  surfaces  and  are  critical  to  explaining 
the  appearance  of  rows  I  and  III  as  lines  of  maxima 
with  2gq  spacing.  The  tetramers,  which  constitute  row  II, 
have  a  more  complicated  geometry  that  gives  rise  to  the 
internal  structure  visible  in  both  the  experimental  and 
simulated  images.  In  the  equilibrium  surface  geometry 
shown  in  Fig.  2,  the  tetramers  are  essentially  symmetric 
(no  buckling).  In  the  course  of  relaxing  from  slightly 
different  initial  geometries,  however,  we  found  a  variation 
on  this  structure  (also  at  equilibrium  and  with  essentially 
the  same  total  energy)  with  a  highly  buckled  tetramer 
(Az  =  0.30  and  0.15  A  for  the  dimer  pair  and  step- 
edge  atom  pair,  respectively).  The  simulated  images  in 
Fig.  3  are  from  a  surface  with  buckled  tetramers;  the 
similar  total  energies  suggest  the  possibility  of  different 
metastable  geometries  within  row  II,  explaining  its  variable 
appearance  in  the  experimental  images. 

Experimentally,  we  observe  that  Si  wafers  oriented 
within  0.5°  of  (114)  form  large  terraces  of  (1 14)  separated 
by  steps.  Wafers  with  orientations  a  few  degrees  off 
(114)  toward  (111)  form  grooved  surfaces  consisting  of 
(114)  and  (113)  facets  [17],  consistent  with  earlier  x-ray 
scattering  results  [5]  and  ex  situ  atomic  force  microscopy 
[5].  These  observations  suggest  that  the  surface  energy  of 
(114)  is  quite  low,  perhaps  comparable  to  Si(l  13)-(3  X  2). 
We  have  investigated  the  stability  of  our  Si(114)-(2  X  1) 
reconstruction  model,  relative  to  low-index  reconstructed 
surfaces,  by  calculating  surface  energies,  Es  =  [ E,  ( /V )  — 
(V£,bulk]/2.  Here  E,(N )  is  the  total  energy  of  the  supercell 
containing  N  atoms,  £,fbulk  is  the  total  energy  per  atom  of 
bulk  Si  (computed  using  the  same  lattice  vectors,  energy 
cutoff,  and  zone  sampling  as  for  the  surface  calculations), 
and  the  factor  of  2  accounts  for  the  two  surfaces  per  unit 
cell.  We  find  the  surface  energies  for  the  fully  relaxed 
2  X  1  reconstructions  of  Si  (111),  (001),  and  (114)  to 
be  84,  86,  and  85  meV/A2,  respectively,  indicating  that 
Si(114)-(2  X  1)  is  approximately  as  stable  as  the  low- 
index  (11 1)-(2  X  1)  and  (001)-(2  X  1)  surfaces. 

While  our  proposed  model — a  terraced  surface  with 
Db  steps  alternately  rebonded  and  nonrebonded — was 
guided  by  our  STM  results,  we  have  also  carried  out 
fully  relaxed  total-energy  calculations  for  a  number  of 
alternative  structures.  All  lead  to  simulated  STM  images 
in  strong  disagreement  with  experiment,  and  most  lead 
to  significantly  higher  surface  energies  than  our  proposed 
model.  Two  of  these  models,  however,  shed  light  on  the 
mechanism  that  ultimately  stabilizes  Si(  1 14):  (model  A)  a 
terraced  surface  with  every  DB  step  rebonded;  and  (model 
B)  a  terraced  surface  with  every  DB  step  nonrebonded. 
Because  rebonded  steps  have  fewer  surface  dangling  bonds 
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than  nonrebonded  steps,  one  might  naively  predict  model 
A  to  be  strongly  favored  relative  to  both  model  B  and 
our  model.  When  fully  relaxed,  however,  we  find  that 
model  A  has  a  surface  energy  of  87  meV /A2,  very  slightly 
higher  than  our  model  (model  B  is  significantly  higher,  at 
97  meV/A2). 

We  believe  that  the  alternation  of  rebonded  and  nonre¬ 
bonded  step  edges  is  ultimately  favored  because  it  leads 
to  an  optimal  balance  between  surface  dangling  bond  re¬ 
duction  and  surface  stress  relief.  For  isolated  DB  steps  on 
Si(001),  rebonded  steps  are  energetically  preferred  to  non¬ 
rebonded  steps  because  the  number  of  dangling  bonds  is 
reduced.  However,  rebonded  DB  steps  contribute  a  ten¬ 
sile  stress  in  the  direction  perpendicular  to  the  step  edge. 
For  sufficiently  wide  terraces,  this  stress  is  elastically  re¬ 
lieved,  so  that  rebonded  steps  are  favored.  For  very  short 
terraces,  however,  this  relief  mechanism  is  not  available, 
and  the  resulting  energy  penalty  may  be  larger  than  the  en¬ 
ergy  gain  from  dangling-bond  reduction.  Si(114)  balances 
these  competing  effects  by  incorporating  both  a  rebonded 
and  a  nonrebonded  DB  step  in  each  unit  cell,  allowing  for 
both  dangling-bond  reduction  and  stress  relief  within  the 
unit  cell. 

This  mechanism  of  stress  relief  can  be  demonstrated 
quantitatively  by  calculating  the  LDA  surface  stress  tensor 
for  each  model,  cr(/-  =  A]  dEj  de,j,  where  e,;  is  the 
surface  strain  tensor  and  A  is  the  surface  cell  area  [18]. 
Based  on  the  above  arguments,  for  model  A  one  expects 
a  large  tensile  stress  in  the  x  direction  (i.e.,  perpendicular 
to  the  step  edge),  whereas  for  model  B  and  our  model  one 
expects  a  smaller  (compressive  or  tensile)  stress.  Indeed, 
we  find  the  values  of  crxx  to  be  0.27  eV/A2  (tensile)  for 
model  A,  —0.03  eV/A2  (compressive)  for  model  B,  and 
0.11  eV/A2  for  our  model.  These  numerical  results  are 
consistent  with  the  qualitative  scenario  described  above: 
model  A  has  a  low  density  of  dangling  bonds  but  high 
surface  stress,  which  incurs  a  large  energy  penalty;  model 
B  has  a  low  surface  stress  but  high  density  of  dangling 
bonds,  which  also  incurs  a  large  energy  penalty.  Our 
proposed  model  compromises  on  both  dangling  bonds  and 
surface  stress,  and  by  so  doing  is  the  most  favorable. 

The  competition  between  the  energy  gain  from  dangling 
bond  reduction  and  the  energy  penalty  from  tensile  stress  at 
rebonded  DB  steps  provides  a  framework  for  understand¬ 
ing  the  surface  morphology  of  a  surprisingly  large  class 
of  orientations,  from  slightly  vicinal  Si(001)  up  to  Si(  1 14) 
( 6  =  19.5°).  For  orientations  up  to  (114),  the  stress  at  re¬ 
bonded  Db  steps  can  be  relieved  either  on  the  terraces  or 
by  the  presence  of  nonrebonded  steps,  and  consequently 
the  terrace-plus-step  morphology  is  favored.  For  exam¬ 
ple,  curved  surfaces  with  local  (119)  and  (117)  orientations 
(6  =  8.9°  and  11.4°,  respectively)  form  well-ordered  B- 
type  terraces  separated  exclusively  by  rebonded  DB  step 
edges  [11,19].  Moreover,  planar  surfaces  with  (115)  ori¬ 
entation  (0  =  15.8°)  consist  of  a  complicated  arrangement 
of  B-type  terraces  separated  by  both  rebonded  and  nonre¬ 
bonded  steps  [17].  Beyond  (114),  however,  the  terraces 


become  too  narrow  to  afford  sufficient  stress  relief,  and 
so  the  terrace-plus-step  morphology  is  no  longer  favored, 
leading  to  fundamentally  different  surface  morphologies. 
For  example,  surfaces  with  orientations  between  (114)  and 
(1 13)  form  mesoscale  sawtoothlike  grooves  [5,9],  and  sur¬ 
faces  with  (113)  orientation  are  stabilized  by  a  more  com¬ 
plex  3X2  reconstruction  [6,7].  We  therefore  conclude 
that  Si(  114)  marks  the  end  point  of  a  family  of  orientations, 
each  member  of  which  has  a  related  surface  morphology 
consisting  of  (001)  terraces  separated  by  DB  steps. 
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